Symmetry breaking in small rotating cloud of trapped ultracold Bose atoms 
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We study the signatures of rotational and phase symmetry breaking in small rotating clouds of 
trapped ultracold Bose atoms by looking at rigorously defined condensate wave function. Rotational 
symmetry breaking occurs in narrow frequency windows, where the ground state of the system has 
degenerated with respect to the total angular momentum, and it leads to a complex wave function 
that exhibits vortices clearly seen as holes in the density, as well as characteristic vorticity. Phase 
symmetry (or gauge symmetry) breaking, on the other hand, is clearly manifested in the interference 
of two independent rotating clouds. 
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Symmetry breaking in finite systems has been a sub- 
ject of intensive debate in physics, in general (cf. the 
Ref. 1]), and in physics of ultracold gases in particular 
over the years. For Bose-Einstein condensates (BEC) 
two symmetries play a particular role: (7(1) phase sym- 
metry and SU(2) (or 5*0(3)) rotational symmetry. In 
the large N limit, one breaks these symmetries by hand, 
as proposed originally by Bogoliubov @- Thus, the accu- 
rate way to deal with macroscopic Bose Einstein conden- 
sates (BEC's) is by the use of a classical field, also called 
an order parameter, or the wave function (WF) of the 
condensate. This function is a single particle (SP) wave 
function, which is the solution of the Gross Pitaevskii 
(GP) equation within the mean field approximation, that 
characterizes the system in a proper way Q- It has an 
arbitrary, but fixed phase, and for rotating systems with 
more than one vortex it exhibits arbitrarily places, but 
fixed vortex array. For dilute ultracold Bose gases (i.e. 
when n|a| 3 << 1 Q where n is the density and a is the 
s-wave scattering length) mean field, or Bogoliubov ap- 
proach is capable to reproduce very well the main proper- 
ties, despite the fact for finite, fixed N and total angular 
momentum L, which are both constants of the motions, 
mean field theory cannot be exact. This observation has 
stimulated a lot of discussion about the nature of the 
phase of BEC 0, , and particle-conserving Bogoliubov 
approach 0. The modern point of view (for a recent 
discussion see 0) implies that two BEC's with fixed N 
each one, will produce a well defined interference pattern 
of fringes as a result of the measurement in only one shot 
(comparable with the calculated n-correlation function) 
in contrast with the density, which would be obtained 
as a mean image of random interference patterns from 
several shots. The position of fringes in the given mea- 
surement are determined by subsequent localization of 
atoms arriving at detectors; the first atom is completely 
random, second is correlated, third even more correlated 
etc. 0. Thus the information about the pattern is ob- 
tained from the many-body wave function by looking at 
pair, triple, ... correlations. The breaking of rotational 
symmetry should occur in large rotating clouds in the 



similar way, and a pure L-state would show, in a time- 
of- flight experiment, a definite interference pattern ac- 
curately represented by n-correlation functions, different 
from a circular symmetric profile of the single particle 
density. It would be a test of the meaning assigned to 
the measurement. Unfortunately, for large N-systems, 
the total angular momentum of the stationary states is 
not well defined and there is no qualitative difference be- 
tween density and n-correlation function, usually showing 
in both cases vortex arrays. For small rotating clouds the 
situation is, however, different, as we have shown in Ref. 
[3- Typically, the GS's are pure-L states for most of the 
values of f2. Only, in the very narrow window of frequen- 
cies, where the ground states is degenerated with respect 
to L, vortex arrays can be obtained, arbitrary small sym- 
metry breaking deformation of the trap potential leads to 
the appearance of symmetry breaking vortex arrays both 
in density and pair-correlations. Namely, in the regime of 
pure L-GS small systems would provide a suitable test for 
the meaning of the measurement distinguishing between 
the density or the pair-correlation output. 

In this Letter we study the effects of symmetry break- 
ing in small rotating clouds of trapped ultracold Bose 
atoms in more depth, by introducing the rigorous defini- 
tion of the condensate wave function, defined as an eigen- 
vector of the one body density matrix operator (OBDM), 
corresponding to the largest eigenvalue. Such definition 
of the order parameter has been introduced in classical 
papers on off-diagonal long range order [T^j. It has, how- 
ever, rarely been used since its application requires the 
knowledge of the full many-body wave function, or at 
least of the exact OBDM. Since for quantum gases exact 
analytic solutions are either not known (2D and 3D), or 
very difficult to handle (ID), so far this definition has 
been only applied to the case of model system with har- 
monic forces. Here we apply for the first time to the 
rotating gas, using exact numerically calculated OBDM 
for few atom systems. We identify in this way possi- 
ble states with vortices, and obtain phase characteris- 
tics of the wave function (reflecting quantized circulation 
of vortices), and provide unambiguous definition of the 
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degree of condensation. With such calculated order pa- 
rameter we then reproduce the density and interference 
patterns for two condensed clouds, and shed new light 
on the discussion of the origins of symmetry breaking in 
finite mesoscopic systems. 

We consider a two-dimensional system of few Bosc 
atoms trapped in a parabolic rotating trap around the 
z-axis. The rotating frequency Q, is strong enough to 
consider the Lowest Landau Level regime with atoms in- 
teracting via contact forces. Our main goal is the descrip- 
tion of the stationary states for different values of f2, an- 
alyzed from the rotating frame of reference, unless other- 
wise stated. Our analysis is performed using the exact di- 
agonalization formalism, valid for arbitrary interactions 
and densities. However, in contrast to the mean field ap- 
proach, this method deals with multi-particle WF's and 
loses the intuitive picture provided by the mean field or- 
der parameter. Our goal is to obtain in a precise way a 
complex scalar field that models efficiently the system, 
and allows to reproduce the important features, such as 
the vortex states. In the regime of relatively low rota- 
tion frequency, where the degree of condensation is high 
and some vortices appear distributed in an ordered ar- 
rays, this scalar field plays the role of a genuine order 
parameter. On the other hand, it looses its capability to 
represent the system as fi approaches the melting point, 
where the prediction is that the vortex lattice dis- 
appears and the systems turns, for large systems, into a 
Laughlin liquid. 

The way to know if there is a "macro-occupied" SP 
wave function in the ground state \GS) is to look at the 
eigenvalues and eigenvectors of the OBDM 0, ^| . That 
is to say, one must solve the eigenvalue equation 

J dr'n {1) (r,r')^(r') =m^{r), (1) 

where 

n W(r,P) = (GS | &(f)$(?)\GS), (2) 

with '5 being the field operator. If there exist a relevant 
eigenvalue ni >rij for / = 2, 3, . . ., then 

VmMry* 1 (3) 

plays the role of the order parameter of the system, 
where 4>i is an arbitrary constant phase. The WF may 
be expanded in the form ipi(r) = y~!™ Pufli^i us- 
ing the complete set of Fock-Darwin (TJ WF's given by 

(pi(r) — e' d9 r l e~ r / 2 /VttU, where I labels the single par- 
ticle angular momentum, and m is equal to the largest 
total angular momentum L involved in the expansion of 
the GS; length unit is here A = y/h/(mu)±), and ui± 
denotes the trap frequency. The same SP basis is used 
in our numerical simulations to represent both the field 
operator and the multiparticle GS wave function. 

An alternative, and perhaps even more appropriate SP 
basis is determined by the functions ipi(r). One can de- 
fine a set of canonical creation and annihilation operators 



for them: 

a] = J dPip* {?)%{?), (4) 

and di being the hermitian conjugate of aj. The Hilbert 
space attain then a tensor structure with respect to the 
modes di, and the new Fock (occupation number) many 

body basis |m) ® |n 2 ) ® The macro-occupied mode 

contains on average n% atoms, but this number fluctuates, 
and most presumably normally, i.e. the fluctuations of 
ni are of order yjn{\ to this aim one has to calculate 
\GS). This implies that atom number fluc- 
tuations between the macro-occupied mode (condensate) 
and the rest of the modes (that could be regarded as 
phonon modes, quasi-particles) will tend to reduce the 
fluctuations of the phase. A natural consequence of this 
observation is to expect that a very fine approximation 
of the GS is given by the coherent state such that 
dx \a\) = Y^Vie 1 ^ 1 \ a i)- If n l for Z = 2, 3, . .. are very 
small we may neglect them, and approximate the many 
body wave function by \ati) <8> IO2) <8> IO3) £g> . . .. In a more 
precise description, we should rather approximate the GS 
by \ai) <g> \a 2 ) <S> \a 3 ) <g) . . ., where di\m) = ^/n'l^ie 1 ' 1 ' 1 \ai) 
where the phases (f>i are arbitrary; one should, however, 
choose them to be random in order to reproduce (on av- 
erage) the same OBDM as the one obtained by exact 
numerical diagonalization. 

This representation implies that the next simplifiying 
step would be the representation of the GS by a clasi- 
cal field entering the GP equation, and containing all 
the involved coherent states \ak), k = l,...m+ 1 as, 

^(r) = Y^k=i ^/^k'4 > k& l ^ k with random phases. Cal- 
culation of quantum mechanical averages would then in 
principle require averaging over random phases, which 
makes this approach technically difficult. 

As long as the exact GS is a state with well defined an- 
gular momentum, (a pure L-state) not degenerated with 
other lowest energy states in different L-subspaces, it is 
easy to demonstrate that the FD functions are the eigen- 
states of Eq.(l) and the eigenvalues n; are the occupa- 
tions usually used in literature. However, at certain val- 
ues of £1 where degeneracy takes place and vortex states 
without circular symmetry (except the case of only one 
centered vortex) are possible Q, the eigenfunctions of 
Eq.(l) are linear combinations of the FD functions and 
the macro-occupied function ipi that represents the vor- 
tex state has expected SP angular momentum given by 
hi = J2j I @ij 1 2 hlj where lj are integers. 

A convenient definition of the degree of condensation 
which senses the loss of macro-occupation is given by 



where N is the total number of atoms and n is the mean 
occupation calculated without the first value n\. 

In what follows, we show some results that confirm the 
convenience to represent the whole system by ipi at cer- 
tain values of fL As a general result, for vortex states, ri\ 
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is always larger than the occupation of the most impor- 
tant FD state within the exact GS. In addition, tpi pro- 
vides a non-ambiguous way to characterize vortices, not 
only showing dimples in the density profile, but also indi- 
cating the position of each one by the change on multiples 
of 2n of the phase S(r) in ipi(r) —\ ipi(r) \ e tS ^ when 
moving around each vortex. In Fig.l for N = 6 atoms, 
we show for three different values of f2 where degeneracy 
takes place, the comparison between the contour plots 
of the density of the exact GS and the density of ipi , as 
well as the map of the phase S(r) of ipi- In the first case 
(a) the GS contains two vortices that appear in a clearer 
way in the order parameter, as it excludes the non con- 
densed part that smears the structure of the GS. The 
same picture is shown in (b) where four vortices become 
visible. In the second case, the map of the phase not only 
localizes vortices with one unit of quantized circulation, 
but also indicates that incipient vortices are growing at 
the edge of the system. In the last case (c), a six- fold 
symmetry is obtained not attached in this case to vor- 
tices, but to a mixed structure of dimples and bumps, a 
precursor of the Wigner type structure observed for few 
atoms in the Laughlin state at an angular momentum of 
L = 30 The degree of condensation as defined in 
Eq.(5) decreases as 0.343, 0.192 and 0.015 from (a) to 
(c). The order in vortices and disorder in atoms evolves 
to order in atoms. As f2 approaches the frequency of the 
trap, the occupations tend to equalize and in the Laugh- 
lin state, where ni are the FD occupations (since it is a 
pure L-state), and the degree of condensation tends to 
zero. 

Some excited states with large L can also be analyzed. 
For N — 3 and L — 9 we obtain a large vortex state 
with three units of circulation. Such state has been pre- 
dicted in previous theoretical studies as a possible giant 
vortex GS (with all vorticity confined to the center of the 
condensate), in the presence of a small quartic potential 
added to the parabolic trap. In such a case stationary 
states for il > uj±_ are possible ^3|. In our calculations 
the giant vortex appears as an excited state, anticipating 
this possibility. So far there is no experimental evidence 
of giant vortex structures in bosonic systems [l4l. they 
have been reported in superconductive disks |l5j . 

Finally we show the interference pattern produced by 
the overlap of two initially independent condensates rep- 
resented by ipi functions. This study is motivated by an 
increasing amount of recent work revealing the possibil- 
ity of obtaining very detailed experimental information 
on the interference pattern produced not only during the 
overlap of two, or more independent condensates |l6[ . but 
also within a unique condensate |17| . 

The idea underlying our assumption is the following: 
we represent the two independent condensates which we 
call a and b by their macroscopic occupied function ip a 
and ipb respectively. By this we mean that the conden- 
sates are in two unknown coherent states \a a ) and \ab) 
from which we know their order parameter except for 
their constant phases 4> a and 4>b (see Eq.(3)). At time 
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FIG. 1: For N = 6 the first two pictures on each row show the 
density contour plot of the GS (p(x,y)) and the ?/>i function 
(pi(x,y)) respectively. The third picture shows tha map of 
the phase S(r) (see text), (a) shows a two vortex structure 
at Q = 0.941 (where degeneracy between L = 10 and 12 
takes place), (b) shows a four vortex structure, £1 = 0.979 
(degeneracy between L = 20, 22 and 24). (c) shows a six-fold 
structure, £1 — 0.983 (degeneracy betwee L = 24, 26, 28 and 
30). In all cases uj± = g = 1 in units of A and u — huj±. 



t = s the condensates are separated by a distance 
d and the traps are switched off. The time evolution 
of the system is obtained (once the transformation to 
the laboratory frame of reference is performed, multiply- 
ing the functions by exp(—iQtL z )) in three steps: First, 
the Fourier transform of the total order parameter (the 
sum of the two contributions) is performed. Then, the 
time evolution of the Fourier components by multiply- 
ing them by exponentials of the type exp(ihk 2 t/2m) is 
realised; this step is done under the assumption that dur- 
ing the time-of- flight the interactions are irrelevant. Fi- 
nally, in the third step, inverse Fourier transformation is 
performed. The results are shown in Fig. 2 where three 
different times are considered. Fortunately, the uncertain 
about the phase relation cf> — <f> a — (f>\> is not important in 
the case considered, as only two terms are involved and a 
change on the relative phase would only produce a global 
shift of the interference pattern. 

We conclude that, we have demonstrated that the use 
of the eigenfunctions of the OBDM operator provides 
a useful and precise tool to analyze the exact GS ob- 
tained from exact diagonalization and specially the vor- 
tex states. These eigenfunctions localize and quantize the 
vortices and reproduce the time evolution of the interfer- 
ence pattern of two overlapping condensates. We want to 
point out that our results imply an alternative interpre- 
tation about a subject that has attracted much attention 
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FIG. 2: Time evolution of the interference pattern during 
the overlap of two released condensates initially separated by 
a distance d = 15A. Initially each condensate contains N = 6 
atoms and their GS are characterized by L — 6 at Q, = 0.019 
and by a mixture of L — 6, 8 and 10 at Q — 0.0847 respectively 
(all quantities are in units of A and u). 



recently related with the interference pattern formation. 
One possibility suggested by Mullin and collaborators @ 
is that the experimental measurement projects the initial 
condensates in Fock states into phase states, the atom 
distribution between the two components become uncer- 
tain and the pattern formation is possible. The other 
possibility discussed by Cederbaum et al. 01 > is that 
the interference pattern appears if one includes interac- 
tion during the time-of-flight even for states that initially 
are Fock states. In our case, the real initial states are 
Fock states and no interaction is included during the 
time-of-flight. However, we assume that the degree of 
condensation of the initial states is large enough to be 
properly represented by an order parameter (condensate 
wave function) . Fluctuations of the number of condensed 
atoms reduce the phase fluctuations and determine the 
order parameter phase. In effect, exact ground state 
manifest themselves as phase states even for small num- 
ber of particles, and in this way the interference patter is 
produced. Note, however, that in our picture the process 
of determination of phase is itself random, and various 
phases <fik are expected to show up from shot to shot. 
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